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Ie INTRODUCTION 


In the Aeronautical Engineering curriculum graduate 
ievel aerodynamics course, AE-4501, the students are exposed 
+o two computer scrograms. One of thes2, prepared by the 
Douglas Aircraft Company, analyzes the potential flow arcund 
three-dimensicnal wings but is ilimited to incompressible 
flow fRef. 1}. The other progran, prepared by Cebeci, 
calculates the friction drag for two dimensional 
incompressible flow over airfoils [Ref. 2]. A serious 
defect of these pregrams is that they are not 
state-of-the-art computer programs. The Douglas program 
dces not consider th2 effects of compressibility and the 
bourdary layer program, in additicn to being restricted to 
incompressible flow, does not predict the laminar to 


muGbulen: transition point. 


A. EACKGRCUND 

In 1980 the Department of Aeronautics azn the Naval 
Postgraduate Schcol acquired an intricat2 computsr program 
recently develoved by the Boeing Commercial Airplane 
Comoany. This state-of-the-art program calculates 


three-dimensional transonic flow over wings and bodies in 


both the cuter-inviscid flow regicn governed by the 
transonic vcoteéential equation and the thin layer in which the 
merst order, cempressible boundary layer equations aré 
assumed +c Le valid. 

The Boeing yprodqram as received was designed to be 
executed on a CDC 6600 or a CYBER 175 computer and was 


written using CDC FORTRAN IV extended language. This thesis 


cr 


herefore was primarily cencerned with th? conversion of the 
orogram to FORTRAN IV extended cempatible with the Naval 
Postgraduate School's (NPS) IBM 3033 systen. The large 
modular voregram was divided so that th potential flow 
analysis ~vcrticn could be run separately. Simplified 


instructions for use of the program were also prepared. 


ie VISCOUS/INVISCID SYSTEF OF PROGRAMS 

The Viscous/Inviscid Wing System (VIWS) of programs 
calculates three-dimensional transonic flow over wings and 
wing kody ccmbinations including details £ «he laminar or 
BeeoUlent flow in the threé-dinensional viscous boundary 
layer. The flew field is calculated in two overlapping 
regiens: an cuter inviscid flow region governed by the 
transonic fotential equation, anda thin boundary layer in 


which the first order, three-dimensional, compressible 
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boundary laver egquaticns are assumed to hold and in which 


he effects of surface téat and mass transfer can be 
computed. A list of the VIWS of programs is presented in 


Table I. 


TABLE I 


Viscous/Inviscid Wing Syst2m of Programs 


FLGz? Jameson-Caughey inviscid, 
+ransonic wing code 


AW11IN Reads gecmetry & velocity. 
data, constructs coordinate 
systen 


VWIN Petential ficw boundary 
layer interface 

AYT1AC1 Three-dimensioral boundary 
layer pregram 

INTERP Peace layer potential 
flow interface 


4Y 
4 
4 


aonb oo) aad 
ah oh 
td'dtn 
(AK) 


Graphics display programs 


The basic sequence of calculations used by the VIWS to 
Obtain matched viscous and inviscid solutions consists of an 
iterative loot in which tke inviscid outer flow analysis and 
she beundary layer analysis are performed sequentially. The 
iterative sequence is ccntinued until either convergence 
(satisfactcry matching) is achieved, or *he naximum numbsr 
of iterations spéecified by the user has besn performed. The 
VIWS programing sequence is shown schematically in Fig. 1.1. 
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Final boundary 
layer calc? 


Generate boundary layer grid 
Ist interface program 


Calculate final boundary layer flow 
1, Boundary layer program—upper surf. 
2. Boundary layer program—lower surf, 


Calculate transonic potential flow 
FEO 27 
Generate boundary layer grid 
Ist interface program 


Calculate 3-D boundary layer flow 
1. Boundary layer program—upper surf. 
2. Boundary layer program—lower surf. 


Plotting programs 


Set up 5° for next iteration 
2nd interface program 

1. Interpolation 

2. 6° under-relaxation 
3. Convergence test 


6° 


converged 


Yes 


Figure 1.1. Viscous/Inviscid Interaction Procedure 


The potential flew is calculated for the bare wing during 
the first iteration. In subsequent iterations, the effect 


of the boundary layer flcw on the outer inviscid flow is 
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felt as a modification to the wing shape through czhe 
addition of the boundary layer displacemant thickness. 
Convergence is recognized, and the iterations are stopped, 
when the maximum change between the new and old displacement 
thickness, expressed as a fracticn of the maximum 
displacement thickness, is less than the convergence 
*olerance chcsen by the user. 

The VIWS utilizes the Jameson-Caughey transonic inviscid 
mero program FLOZ/, ¢t0O carry cut the potential flow 
analysis. The boundary layer analysis is performed by a 
finite difference bcoundary layer prediction program 
d2veloved by the Being Cecmmerciai Airplane Company. The 
basic ztheorv Fehind the bcundary layer program is contained 
in [Ref. 3]. A detailed descripticn of th2 VIWS of programs 
(excluding the potential flow program FLO27) is contained 
in (Ref. 41. A basic guide to the use of the VIWS of 


programs is contained in [{Ref. 5]. 
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TI. POTENTIAL FLOW PROGRAM FLO27 


Eecause of the extensive length and number of program 
modules in the VIWS, the Fotential Flow Program, FLO27, was 
singled out for cenversion. It was anticipated that FLO27 
would be run separately at first and racombined with the 
other croqram modules at scme later date when these modules 
were themselves translated for execution on the IBM 3033 


computer. 


A. RE~PROGRAMING 

The Potential Flow Preqran, hereafter called FLO27, was 
received cn magnetic tepe and loaded into «he IBM 3033 mass 
stcrage system using the Jcb Control Languéeg2 (JCL) routines 
presented in Arpendix A. The magnétic tape contained twenty 
(20) total files in which «he format was 9 track, 1600 CPI, 
unlakled. The card image format for the sixteen (16) 
program files is 80 characters per record and the four (4) 
Output files ccntain 150 characters per record. The program 
and cutout files on the criginal CDC tape ars listed in 
faole Il. 

The FLO2Z?7 pregram was converted to FORTRAN IV extended 


Suitable for execution on the IBM computer using the NPS CDC 
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TABLE II 


CDC Magnetic Tape Filés 


lie az 3G ELOi27 Potential Flow Program 


z /3194 AYTIIN Reads gecmetry & velocity data 
constructs coordinate systen 

s47 378 VWIN Potential flow boundary-layer 
interface 

4 /6611 A411A01 ee ce dimensional boundary 
lay#r program 

SY oo INTERP FOES CSI hoe potential flew 
interface 


3770 | cours 

: 

Procedures files 
Boeing McLean computer progran 
Pie _| ourire [output from VIN 


OUT4YII1L | Output from boundary-layer, 
lower surface 
OUT411U | Output fErom bdoundary-layer, 


upper surface 
<O ITEM conversicn guide [Bef. 6]. The first step taken 


consisted cf crogram compilation uSing the WATFIV compiler 


with its extended error messages. The listing which was 
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produced flagged all areas of the program which required 
révision. Pregram changes were accomplished utilizing *his 
WATFIV listing. Seme of the mcre general and repetitive 


changes are listed in Table III. 


TABEE OL Yi 


FLO27 Re-Programing Changes 


cDC Code IBM Code Change 


Variables: FREAD, FREAF, a ! 

Peel t, EWRIF, TREAD, Eliminated from program 
ITREAF, IWRIT, IWRIFP- 

WRITE ( IWRIT,600) WRITE (6 ,600) 
READ(IREAF,500) READ (5,500) 


READ ee WRITE 7 or Changed to READ 14, 
REWIND 7 WRITE 14 or REWIND 14 


All of this type esliminated 


The most difficult change +o make occurred with the CDC 
Buffer IN or Buffer OUT statements which were used in the 
program to transfer portions of a three-dimensional array 
into and out cf main memory. The Buffer routines reduce tkh2 
memory size required to execute the program. This statement 
tvyce cccurred inthe main program and several of the 
subrcutines. 
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The change reguired to translate this statement is presented 
below with the CDC code preceeding the IBM FORTRAN. 


EUFFER OUT (N3,1) (G(1,1,1),G(MX,MY,1)) changed to 


WRITE (N3) ((G(I,J,1) -I=1,MX) ,J=1,MY) and 
POreeR IN (N1,1) (6(1,1,8),6 (UX, MYM) } Charged to 
READ(N1,ERR= ) ((G(1,J,4) ,1=1,MX) ,J=1,MY) 


Tht variable ERR was assigned the GO TO Statement number 
of the UNIT statemen* immediately following the BUFFER IN 
Mere of code. AS an example, if the UNIT statement 
follcwino the EUFFER IN cede was - If (UNIT (N1).GT.0.) GO TO 
151, ‘then the number 151 was assigned to variable ERR 
follcwing the equal Sign. All CDC UNIT statements were 
eliminated from thea FLO27 source code per Table II. 

In addition to the precgram changes required to run FLO27 
on the IBM ccmcuzer, several lines of codes were added to 
modify the outvut format +o a more usable form. A 
subrcutine, VERTEC, which calls the Versatec plotter was 
also maacd to enhance pregram usefulness. Thus plocting 
routine is user controlled through an input variable and is 
ex lained in the next section. The modified FLO27 prcgram 


source code is presented in Appendix E. 
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To facilitats program data entry several input variables 
which had recemmended values were initialized to these 
values within the Main program and the subroutine GEOM. The 
initialized inpu* variables and their values are presented 


in Table IV. 


TABLE IV 


Initialized Input Variables 


[scan | 
2 
[recone fo 
TC 
sice 
[—erce 
sme 
a a 
eT 
a a a 
[_astus fit 
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A comelete description of each input variable in Table 


IV can be found on pages 19 through 23 of [Ref. 5]j. 


B. PROGRAM DESCRIPTION 

The FLC27 program is a computer codé written +o analyze 
+he transcnic flow over awing alcne or a wing ona 
cylindrical fuselage. It uses a finite-volume formulation 
=o sclve the exact potential flow equation in conservative 
OCT « In the dévelopment of the equations, the basic 
aSSumftions are; sore flow, no heat or work transfer, 


isentropic flow, irrotational flow, a2o body forces and a 


verfect gas. The velocity vector in cartesian coordinates 
is 
aN “> “N 
V= ul + vi + wk C2) 
whare u, vand w are +tfe velocity components. The 


continuitv equation, assuming steady flow, is 


os (eu) eter ew) = 0 (22) 
Next a velocity potential is introduced such that the 


velocity components are calculatéd as the gradienz of this 


notential. 


u =9 ,v =O 7 w =@ (2.3) 
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With the introduction of the velocity potential, the 


continuity equation 2.2 becomes 
2h) + 0H) + Shem) = 0 (2.4) 


Assuming no heat or work transfer, the @nergy equation can 


be written as 
2 
os (255) 


The flow is assumed to be uniform in the far field. On the 
surface of the bedy, the normal velocity component is zero. 
The velocities and densities of the near field are 
normalized using the free stream velocity and density, <hus 
Veo = 1 and boo = t Using the assumptions that the flow is 
isentropic and a perfect gas, the energy equation 2.5 can bea 


shown to -e 


=|, Z -1 
0 = 1+ iy Mo. (-¥) Y (2.6) 


With equaticns 2.5 and 2.6 there ar2 only two unknowns, 
D and P ~ They can be solved, subject tc «he boundary 
condition of flcw tangency, using a finite volume technique. 


Tee basic numerical scheme for <zhe solution is the 
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constructicn cf a mesh frem small volume elements (cubes) 
me 

which are packed around tke wing or wing body configuration. 

The cubes in the computational domain are separately mapped 


-9 distorted cukées in the physical domain by independent 


transformations from local coordinates X, Y and Z to 
Cartesian coordinates x, yandz. The mesh points are the 
yver-ices (ccrners) of the mapped cubes. The velocity 


potential and dersity are calculated at each vertex in the 


mesh. The pressure distribution can then be calculated fron 


= (2.7) 


In the event that tke local flow velocity becomes 
Sucersonic and shocks occur, these are handled in the usual 
manner by insuring that: 

1) The tangential’ velocity components are egual on each 
side of the shock. 

2) Continuity is maintained by keeping the product of 
YU, constant across the shock (where J. is the normal 
velocitv component). 

3) Discontinuous expansions (corresponding to an 


Nexpansicn shock") are excluded from the flow field. 
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The assumoticn cf isentrofic flow along with the existence 
eof shecks presents a contradicticn which can only be 
resolved by limiting the flow to very weak shocks for which 
entrepy and verticity generation may be ignored. Thus, 
solutions will be valid only for subsonic free stream 
velocities. 

The main three-dimensional array containing «the 
votential function data is stored on disk, and special 
unformated incutyoutput statements are used to bring planes 
or data into céntral computer memory and to store updated 
planes of data back on the disk. In th? construction of the 
computaticnal coordinate system, a Joukowski transformation 
is used to transform *he cylindrical fuselage to 2 vertical 
slit and then a sheared parabolic transformation is used in 
planes containirg the airfoil seéctions. A detailed 
mathematical formulation of the potential flow analysis is 
contained in [ Ref. 7]. 


ie ercaram input 


The inrut to FLO27 consists of variables which are 


read with an &F10.6 FORMAT. Each input card has a titl 


iD 


Card which precedes it. ThIis ti*#le card contains the inpu~ 


variakle name and effectively labels the input data for easy 
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reference. The title for each input data is placed in the 
same column as the input data it labels. The +itle cards 
are read with a 20A4 FORMAT. All numerical input valu2s are 
real numbers. The following data deck, listed card by card, 
is the mininum input data required for a simplé¢ wing 
analysis. Each "card" can be interpreted as one line of 
data cn your terminal. A complete sample data set is 
presented in Aoreéendix C. 
CARD 1 The Run Title (64 characters maximum) 
CARD 2 Title card for the input variables 
PNK, ENY, FNZ, FMESH Aha FPLOT 
CARD 3 
Geils 1-10 FNX - Number cf computational cells in 
the chordwise direction for the initial mesh. 
MAX = 160/2**n, wkere n = FMESH -1. (See Cols. 31-40 
fom etic i) 
Cols. 11-20 FNY - Number of computational cells in 
the nermal direction frem the airfoil surface for 
the initial m2sh. 
MAX = 16/2**n, where n = FMESH -1. 
Cols. 21-30 FNZ - Number of computational calls in 
the ssanwise direction for the initial mesh. 


MAX = 32/2**n, where n = FMESH -1. 


Z3 


Cols. 21-40 FMESH - Determines the number of times 3a 
program generated computational mesh is refined. 
EN ceeacwiyet.oe 2.0 OF 320 for coarse, medium or 
fine mesh. If 3.0 is selected the prcegram will 
calculate flow over the wing for the coarse mesh 
then half the mesh size (medium), recalculate, then 
half the mesh again (fine) and doa final potential 
flow calculation. Output parameters are printed for 
each mesh size fer which calculations were 
performed. 

Cols. 41-50 FFLOT - Output flag 


0.0 NOotmal cutput without printez-plot of Cp 


1.0 = Normal output with printer-plot of Co 
at each computaticnal mesh point for 
@ach wing section. 
2.0 = Normal cutput with Versatec plots of 
Cp versus X/C for each wing section of 
the final mesh. 
CARD 4 Title card for the input variables 
FIT, COVO and P10 
CARD 5-M Cne card for each computational mesh. Total 


numker of cards equal to ™ = FMESH. 
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Gomis. 1-10 FIT - A parameter which fixes the 
Maximum number of iterations the program will us¢ to 
cenverge the velccity petential toa specified 
tolerance (COVO). This parameter must be repeated 
for each mesh refinement. 
Cols. 11-20 COVO - Velocity potential convergencs 
Grwtscria. This input variable is also entered for 
each selected mesh. A value of 0.000001 is 
reccmmended. 
Cols. 21-30 P10 - This parameter determines the 
Suksonic point relaxation factor for the specified 
mesh size. A value of less than 2.9 must be entered 
for each designated mesh. Pecommeandsd values are: 
1.6 fer coarse, 1.3 for medium and 1.2 for the fins 
mesh. 

CARD 6 Title card for the input variables 


FMACH, YA, AL and CDO 


Cols. 1-10 FMACH - Free stream Mach number 


Cols. 11-20 YA - Yaw angle in deaqrses 


Cols. 21-30 AL - Angle of attack in degrees 
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Sand &g 


CARD 9 


Genus. 31-80 CDO  - Drag ceetfticisnat due to |ekin 


fee CGE ven > Unless known, an estimated value of 0.01 
is reccmmended. 
Title card for the input variables 


ZSYM, FNS, SWEEP, LIHED and FUS 


Cols. 1-10 ZSYM ~- The wing planform symmetry 


SE lgqer. 
0.0 = Yawed wing, has no sSpanwis? symmetry 
1.0 = Swept wing, has Srfanwise symmetry 
Cols. 11-20 FNS -This input variable <zells the 


cregram the total number of Wing sections you have 

selected to define the wing half span. The number 

must ke at least three (3) but not more than eleven 

(11) sections. 

Cols. 21-30 SWEEF - Leading edge sweep angle in 

deqréees. 

Cols. 31-40 DIHED - Dihedral angl¢ in degrees. See 
Pig «2a 1 

Cols. 41-50 FUS - Input the fuselage radius. Enter 


O.9 for a wing-alcne case, 
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Figure 2.1. Dihedral Angle 


Data input cards from 10 through 15 are used for 
d2fining wirg planforms and section géometrics. For the 
Fest Sea sécticn, all data cards from card 10 through card 
15 must be used. For the second and subsequent sections 
there is an cption for skipping the wing section defining 
data (cards 12 through 15) and copying the data from that of 
the pravicus section. This option is controlled by the input 
Variable FSEC. If this cption is not used, data scards from 
10 through 15 must be repeated for each wing d2fining 
section. The number of wing sections which are defined is 
insut with the variable FNS. Remember, up to 11 sections may 
be defined, anda minimum of 3 sections must be defined. 
All wing clanform and section defining geometrics must be in 


consistent units. Wing planform and section defining 


.wantities are cresented in Fid. 2.2. 
oT 


YLe— chord 4 


1 
Y FN 
X 


Figure 2.2. Wing Defining Geometry 


CARD 10 Title card 


ZS, At, YL, 


CARD 11 


Cois« 


Celis. 


Cols. 


COS. 


d= 20 


21-30 


31-40 


for the input variables 


CHORD, THICK, AT and FSEC 


ZS - The section spanwise coordinates 
(Start at th? canterline and werk 
outboard) 

XL - Section leading edge X coordinets 

YL - Section leading edg2 Y coordinaits 


CHORD - Section chord length 
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CARD 12 


Cels. 41-50 THICK - The thickness scaling factor 


can be used to scale all Y coordinates of tné wing 
section. Thus percent thickness and camber are 
increased (or decreased) eccordingly. Use 1.0 if 
no scaling is desired. 

Cols. 51-60 AT - Th2 twist angle of each section 
(gécmetric twist) measured from the X axis to the 
ehovd = ine . A positive twist angle reduces the 
section angle of attack and gives "washout". Use 
O20 fcr no twist. 

Cols. 61-70 PSEC - This is @ flag which determines 
whether or not the program reads wing section 
definino geometry from a previous wing section or 
from new defining geometry. For the first section 
defined you must set FSEC tc 1.0. Following the 
first sectior, if you define new section geometry 
then use FSEC = 1.0. I£ you want the program +o 
read the section geometry defined from the previous 
s¢ecticn then set FSEC = 0.0. 

Title card for the input variable 


FN 


Zo 


SARD 13 
Celis. 1-10 FN + This variable contains the number 
cf pcints which define the upper and lower surface 
of the section. A maximum of 161 points may be 
used. 

CARD 14 Title card for tke input variables 
XF (I) and YP (I) 

CARDS 15-1 te 15-N Total number of cards ¢guals N, 
where N = integer part of (FN+2)/3. 
The X and Y coordinates at each point are entered 
in pairs, three fpoints to a card. (S2e Appendix C 
fer samvol> input) 
Cas. Wake ey = k& COOTdi netesof the. wing 

section point 
11-20 YP(I) - Y coordinate of the wing 


secticn point 


21-30 defining X coordinate for next 
rPorne 

31-40 defining Y coordinate for next 
point 

41-50 defining X coordinate for 


following point 
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5 1-60 defining Y coordinate for 
following point 

The X and Y coordinates of the wing section defining 
points must be entered starting with the upper surface 
railing edge pecint and preceeding along the upper surface 
to the leading edge, and returning along the lower surface 
to the lower surface trailing edge point. It is very 
im.ortant to define the section leading edge with a lazgs¢ 
number of closely spaced points. Suggest at least 9.05 
snacing or less retween X coordinate values from 0.1 X/C to 
the leading edge, X/C = 0.0. Each X and Y coordinat2 point 
is nermalized using the chord length for that section. 


Secticn defining geometrics are illustrated in Fig. 2.3. 


— chord — 


Figure 2.3. Section Defining Geometry 
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CARD 16 Title card ccntaining the words in Cols. 1-80 


END OF CALCULATICNS 
CARD 17 Title card for the input variable 

FNX 
CARD 18 

Cols. 1-10 FNX - This variable indicates the end 

of a set of calculations and must be set equal +o 

0.0. Its puzpose is to indicat that the program 

has run ts completion. 

2. Pregqram Output 
Outout from the FLC27 program varies with the value 

of the input variable FPLCTI. When FPLOT is set equal to 0.0 
a rormal output is produced. This nermal output contains 
(in crder of occurance): refined input geometry data 
*neluding trailing edge slope andangl2 calculations; 
-=terative soluticn of the potential flow mesh; section 
characteristics and wing characteristics. The iterative 
eS secticn and wing characteristic data are repeated 
fOr each mesh refinement requested. PHYS, Le one input 
Variable FMESH is set equal to 3, these data are calculated 
and cutput three times. The last data in the normal output 


consists cf the non-dimensionalized chord (X/C) and vressure 
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coefficient (Cp) data at each computational mesh point for 


each wing section calculated during the final mesh. A 
sample of the normal output data is presented in Appendix D 
and rerresents the output data from Appandix C input daté. 

If variable FPLOT is set ¢egqual +9 1.0, czhe output 
data is increased considerably. This output contains +he 
normal output pius a line printer-plot of the-«pressure 
coefficient at each computational mesh point for @ach wing 
section. The line printer-plot is producad for each wing 
section of each mesh refinement. The length of the cutput 
data with FPLOT set equal to 1.0 can approach 6000 records 
d2-ending on the number of mesh refinements requested. 
These vlots are of questionable value and, therefore, an 
alternate plotting program was developed. 

When the variable FPLOT is set equal to 2.0, the 
normal output data is preduced plus a Versatéec plotting 
Subrcutine (VERTEC) is called. The subroutine outputs, vie 
the Versatec crlottar, plots of Cp versus X/C for each wing 
section of the final mesh calculaticns. This reutine is 
Simply putting into plot form the Cp and X/C numerical data 
contained in the normal output. A sample of the Versatec 


plozc is soresented in Fig. 2.4. 
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Figure 2.8. VYVersatec Plot of Cp vs. X/C 
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IIT. INTERACTIVE INPUT PROGRAM FLO27IN 


The input data file required for tha PFLO27 program is 
extensive. Errors in input data FORMAT will cause progran 
errors at execution time. In crder to eliminats these 
errors and reduce the input data workload, a computer 
+arminal interactive program was written. This interactive 
program, Called FLO27IN, is a user-frisndly way of creating 
an input data file fcr the potential flow wing analysis 
program FLO27. The FLO2Z7IN program source code is presented 
in Arrendix F. 

The interactive program, FLO27IN, when executed ask's 
questicns cf the user in crder te construct and write to the 
user's "A™ disk the required FLO27 input data file. The 
follicwing presents the step-by-step procedure for executing 
she interactive program FLO27IN. 

STEP #1---Log on to any IBM 3033 interactiv2 terminal 
with your user number and password. 
STEP #2---Once logged on and in the CMS operation mode 
type: 
CP LIK O247P 191 120 RR then hit ENTER: 
STEP #3---The werd PASSWORD will appear, Type and FNTER 


AERO 
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STEP #4---Tyoe and ENTER 


ACC 120 D 
STEP #5---Type and ENTER 
LOAD FLO2Z7IN (START 

The screen will display the header for the interactive 
orog ram. Answer 2ach question presented. At the end of 
each guestion in parenthesis is the input data variable 
associated with that question and whether the input 
rarameter is a real number (R) or an integer (I). Example: 
==> Fnter the free stream Mach number (FMACH): (R). FMACH 
is the invut data variable for the question. As you proceed 
emecough «he FLC27IN progran, Opportunities to review and 
change input data will be presented. Should it become 
necessary to change your input data after completing the 
FLO2Z7IN srogram, you can simply XEDIT the created data file. 

The FLO2Z7IN program also incerforat2es a library which 
contains the wing-section defining data fora number of 
current wing shaces. A copy of this library is presented in 
me endix B. This feature will be displayed during program 
execution by the use of a menu from which the user can 


select a pre-defined wing section or define his own. 
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Upon ccmpletion of user inputs to the interactive 
uroaztam three additicnal data lines aré automatically 
written to the bcttom of the input file. They are: 

END OF CALCULATION 
FNX 
9.0 

ma additicn, Job Control Language (JCL) cards are 
Wwoitten *o the top and bottcem of the file. All JCL cards 
start with a // format. After FLO27IN has run to completion 
ty-e and enter RELEASE 191 to release the aero disk which 
was linked while executina the PLO2Z7IN program. The created 
data file is written to the user's "A" disk with <filename> 
<file+voe> of FLO27 DATAIN. Additional changes can be made 


Simply by entering th2 XEDIT mode and editing the file. 
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iweb Ov me OATLGH SY STEM EXECUTION 


The votential flow pregram FLO2Z7 can be executed after 
the inout data file has been created. The batch processor 
is reguired fcr FLO27 execution because of the extensive CPU 
time needed to run the progran. While in the XEDIT mode, a 
standard JCB card must be added +o the top cf the FLO27 
DATAIN file orior to submission for job 2xecution. The JOB 
card has the form: 

//jotname JOB (nnnn,prepp) ,tident!,CLASS=J 
jobname = may contain up to 8 alphanumeric characters, 
the first of which must b2 alphabetic. 


Enon = your user number 


Ppppp project number, assigned by professor 

tident' = contains the user's own identification 
information. A maximum of 20 characters may be 
contained within the single quotation marks. 

After adding the JOB card to your data file, you are 
ready to execute the program. Type SUBMIT FLO27 DATAIN and 
press ENTER. Eatch runs are normally net worth waiting for. 
To inguire about the status of the job, enter INQ and the 
job name used on the JOB card or "logoff". If the system is 


busy and the maximum m¢esh size was selected, it may be 


several hours before your job is run. 
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When the jok is run the output will be spooled to «the 


batch printer located next to the VM printer in the main 


ct 


computer building. The ticle at the top of the printout for 
batch jobs is the name entered on the JOB card. WE as 
desired to have the program output data spooled directly to 
the terminal, it will be nécéssary to add one additional JCL 
card to the input data set. This card must be placed 
immediately following the JOB card and has th2 form: 
//*MAIN ORG=NPGVM1.nnnnP 
where nnnn = ycur user numter 

Inserting this card in the input data will cause all 
program output to be spocled to the user's virtual reader 
where it may Ee looked at, printed cr transferred to his "aA" 
disk. To enquire as to whether information is in the reader 
Simply type RDR and hit enter, then follow the instructions 


on the screen. 
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Ve. PROGFAM TEST RESULTS 


The FLO27 program was tested in three stages; {1) during 
«he reprograming phase for conversion completeness, (2) 
after sucessful convérsicn with suitable wing data for 


program accuracy and (3) during an AE-4501 class project. 


A. ACCEPTANCE TEST DATA 

To test and ensure that *he FLO27 program was converted 
to ITEM compatiktle Fortran without error, an acceptance test 
data set was used. The acceptance test input and output 
data was supplied with the criginal CDC program source code. 
After conversion of the FLO27 program to Fortran suitable 
for the NES IEM system, the acceptance test input data were 
Tun and che output results compared to the output generated 
by tne CDC systen. 

Poth outsut data sets were numerically exact when the 
FLO27 proqram was run in double preécisicn on the IBM systen. 
If the oregqram was run ir single precision, the numerical 
output values were exact te the third décimal place. The 
difference in single precision accuracy occurs because the 
CDC system uses a 64 bit word length while the IBM systen 


word length in single precision is only 32 bits. It was 
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decided that the IBM single precision accuracy was 


satisfactcry. 


B. CCMPARISCN WITH OTHER FROGRAMS 

The FLC27 program was also tested for accuracy by using 
the wing olanfcrm and section data from 2 NACA 572 wing. 
The data were run on both the FLO27 program and the Douglas 
potential flow program {Ref. 1}. The data generated by both 
programs was ccmcared to wind tunnel data fcr the NACA 572 
wing (Ref. 8]. The results are presented in Pig. 5.1 as 
plots cf lift ccefficient versus angle-of-attack. The 
results shew that for the NACA 572 wing the FLO27 program 
more accurately predicts the wing lift coefficient than does 


the Cocuglas pregran. 


=. AE-4501 CLASS PROJECT 

The final test phase was conducted by introducing the 
FLO27 program into the AE-4501 course as a class project. 
This was accomplished to determine student problems/comments 
concerning the data input program FLO27IN and to test an 
additional wing shape. The wing chosen for study was that 
of the A-7 airplane. The A-7 wing has a distinct leading 
edge netch at the approximate mid-span. When the planfern 


geometry was run with the notch included the FLO27 progran 
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LIFT COEFFICIENT - CL 


6.0 -40 -20 00 20 40 #60 80 100 120 
ANGLE OF ATTACK (DEG) 


Figure 5.1. Program Calculated and Wind Tunnel Data 


Tan to completion but gave negative values for section and 


total induced drag coefficient. The value for the lift 
coefficient was low for the freestream Mach and 


angle-of-attack used. It was found that if the notch was 
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excluded from the wing gecmetry input data the program 


results were satisfactory both for induced drag and lift 
coefficient. 

From the AE-4501 class experience it was determined «hat 
share wing tclanform discontinuities cannot be handled by «he 
progran. If however, the changes in shape are gradual, such 
as awing glovz, the program cutput appears to be 
Satisfactory. Such was the case with the acceptance test 
case data where the wing geometry was ‘that of the F-8 


supercritical wing which incorporates a wing glove. 
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AEPENDIX A 


icient space on *h2 mass 
Pee pec Olc ch ts 
eS = 


C This JCL routins allocates su ft 
= a 
9*, ,CLASS=A 


C stcrage2 svsten to store the en 


7/IACK JOB (2266,0178) , * PASCHALL- 

//*HAIN ORG=NEGVM1. 3266P 

77 EXEC PGM=TEFBR14 

7/DD1 DD UNIT=3230V, MS VG P=PU BEC DISP= (NEW CATLG) 

We DSN=NSS.S3266. WFELOW. LATA,SPACE=(CYL, ( 16,4, 2)) 

Va 

C This JCL, routine is used to transfer all tape files to 
C a partitioned data set in the mass Bee ees systen 
7/JACK JOB (2266,0178) , "PASCHALL-2759' ,CLASS=J 

//*MAIN ORG=NEGVM1. 3266P 

7/COEY EROC FILEE,MEM= 

//  ¥XEC PGM=IEBGENER 

//SYUSPRINT DD SYSOUT=A 

7/SYSIN DD DUMMY 

//SYSUT1 DD UNIT=340 0-6, VCL=SER=WFLOW, DIS P= (OLD, PASS), 
47 LABEL=(SFILE,BLP,,1N) 

//  DCB= (RECFM=F,BLKSIZE=£0, DEN =3, OPTCD=Q) 

//SUSUT2 BD D1SP= (OLD, KEEE) , 

7/ DSN=MSS. $2266. WELOW. SOURCE (6MEM) 

//  DCB= (RECFM=FB, LRECL=80, BLKS 12 F=6400) 

// PEND 

// EXEC COEY,FILE=1,MEM=FLO27 

Y/ EXEC COPY, FILF=2,MEM=A411IN 

// EXEC COEY, FILE=3, MEM=VWIN 

Yf EXEC COPY, FILE=4,MEM=A411A01 

// EXEC COEY, FILE=5,MEM=INTERE 

// EXEC COPY, FILE=6, MEM=A411PS 

//  &XEC COEY,FILE=7, MEM=AU11P 1 

Jf EXEC COPY, PILE=8, MEM=A411P2 

// EXEC COEY,FILE=9, MEM=COUPLE 

// EXEC COPY, FILE=16, MEM=ITER 

// EXEC COEY,FILE=11, MEM=DATAIN 

// EXEC COPY, FILE=12,MEM=PFINAL 

Jf SBXEC COEY,FILE=13, MEM=BOEB 1 

7/ EXEC COPY, PILE=14. MEM=CONTPLT 

// EXEC COEY,FILE= 15, MEM=CORD PLT 

// EXEC COPY, FILE=16 , MEM=STRE PLT 

//CCEY2 FROC FILE=, MEM=, LRECL=80, BLK=6400 

// EXEC PGM=IEBGENER 

//SYSPRINT DD SYSOUT=A 

//SYSIN DD DUMMY 

//SYSUT1 DD UNIT=3400-6, VOL=SER=WPLOW, DISP=(OLD,PASS), 
// LABEL=(SFILE,BLP,,IN) 

// DCB= (RECFM=F 3B LKSLZE=SLRECL , DEN=3 ,OPTCD=0) 
//SYSUT2 DD DiS B= (OLD, KEEF) DSN=MSS.S 3266.7 FLOW.DATA (SMEM) , 
yy DCB= (RECFN=FB, LRE L=SLRECL, BLKSIZE=SBL&) 

// EXEC COPY2,FPILE=17, LRECL=150, BLK=6000,MEM=OUTF27 
f/f BXEC COFY2,FILE=18, LRECL=150, BLK=6000, MEM=OUTIFC 
// EXEC COPY2,PILE=19, LRECL=150,5LK=6000,MEM=OUT4U11L 
44 EXEC COEFY2,FILE=20, LRECL=150, BLK=6000, MEM=O0 T4110 
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APPENDIX B 


PEEbanieOr ATRECTL SECTION GEOMETRIES 


user input section coordinate data 


fla< 


plate cata 


symmetrical 


supercritical wing (cambered, 


chord) 


NACA 24-30-0 (cambered, 
F-14 wing (cambered, 
A-7 wing 
chord) 
LISSAMAN 7769 Airfoil 


43% 


Enorc) 


NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 
NACA 


00 10 


(symmetrical, 


0010-34 (symmetrical, 
0010-35 (symmetrical, 
0010-64 (symmetrical, 
0010-66 (symmetrical, 


tG6=002 
63-070 
63 A010 
64-010 
644010 
ela) 10) 1/8, 
654010 
66-010 


(Symmetrical, 
(Symmetrical, 
(symmetrical, 
(Symmetrical, 
(symmetrical, 
(Symmetrical, 
(Symmetrical, 


(Symmetrical, 
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(cambered, 


10% thickness at 


10% thickness 
10% thickness 
10% thickness 


Wing (11% thickness at 30% chord) 
12% thickness at 32% 


11% thickness at 


12% thickness at 30% chord) 
9.5% thickness 2t 37% chord) 
(7 deg droop at 20% chord, 7% thickness at 


30% 


30% chord) 
10% thickness at 40% 


at 
elie 


at 


5 0% 
40% 
60% 


Chord) 
chord) 
eines d) 


chord) 


9% thickness at 50% chord) 
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10% 
10% 
10% 
10% 
10% 
10% 


thickness 
thickness 
thickness 
thickness 
thickness 
thickné ss 


ehnickness 


35 % 
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“OK 
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